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A HOT SPOT IN COMA 
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ABSTRACT 

We study the temperature structure of the central part (r < 18'~ 0.7ft^ Mpc) of the Coma cluster of 
galaxies using ASCA data. Two different analysis methods produce results in good agreement with each 
other and reveal the presence of interesting structures in the gas temperature distribution. Globally, the 
average temperature in the center of the cluster is 9.0 ±0.6 keV in good agreement with previous results. 
Superimposed on this, we find a cool area with temperatures of 4-6 keV associated with a filament of 
X-ray emission extending southeast from the cluster center detected by Vikhlinin and coworkers. We 
also find a hot spot with a temperature of around 13 keV displaced north from the central peak of 
emission. The distribution of the gas temperatures and relative specific entropies suggests that the cool 
features are most likely gas stripped from a galaxy group centered on NGC 4874 falling toward the core 
from outside, while the hot spot located "ahead" of this in-falling gas is due to shock heating. Thus 
our results suggest that we are observing Coma during a minor merger with a small group of galaxies 
associated with NGC 4874 shortly before the initial core passage. 

Subject headings: galaxies: clusters: individual (Coma, A1656) — galaxies: ICM — X-rays: galaxies 



1. INTRODUCTION 

The rich, nearby Coma cluster (A1656) has long been 
considered as the archetype of a relaxed cluster system. 
However, with the advent of X-ray imaging satellites (Ein- 
stein and ROSAT), as well as extensive optical studies, it 
has become clear that Coma, like many other clusters (e.g., 
Geller & Beers 1982; Jones & Forman 1983; Bird 1994) ex- 
hibits significant substructure at various linear scales. 

At the largest scales (~ 3 ftjjjj Mpc), R OSA T observa- 
tions have shown a large subcluster ~ 50' to the southwest 
of the center associated with NGC 4839 (Briel et al. 1992; 
White et al. 1993; Burns et al. 1994). Optical studies 
( Mellier et al. 1988, Biviano et al. 1996, hereafter B96; 
Colless & Dunn 1996, hereafter CD96) have found sig- 
nificant enhancements in the density of galaxies associ- 
ated with this subcluster. At intermediate^ i00h^ kpc) 
and smaller scales these authors as well as several more 
(Fichett & Webster 1987; White et al. 1993; Mohr, Fab- 
ricant, & Geller 1993) have found significant substructure 
in the X-rays and complex velocity distributions in the op- 
tical and associated with the two dominant galaxies NGC 
4874 and NGC 4889. 

Several authors (Davis & Mushotzky 1993; White 
et al. 1993; Vikhlinin et al. 1994) have found distinctive en- 
hancements of X-ray emission associated with the galaxy 
concentrations around the two dominant galaxies, NGC 
4889 and NGC 4874, for which they obtained mass esti- 
mates consistent with those from the velocity dispersion 
inside these galaxy groups (Mellier et al. 1988). A more 
detailed wavelet analysis of the ROSAT image (Vikhlinin 
et al. 1997, hereafter V97) has also revealed a gas filament 
extending southeast from the cluster core for ~ 1 h^Q Mpc. 
One explanation for this filament is as a trail of gas left in 



the wake of an in-falling galaxy group. 

Temperature maps of the X-ray emitting gas are a new 
and especially sensitive tool for the detection of dynamic 
activity. Hydrodynamic simulations indicate that merg- 
ers should produce characteristic temperature patterns 
that survive 4-6 times longer than perturbations in the 
gas density (e.g. Schindler & Muller 1993; flicker 1997). 
Recent observations (Henry & Briel 1995; Hcnriksen & 
Markevitch 1996; Markevitch et al. 1996b, 1998; Donnelly 
et al. 1998) have found just such characteristic tempera- 
ture structures in a variety of clusters, which have been 
interpreted as due to cluster mergers. 

A large-scale (r = 60' ~ 2.4 /ijjj, 1 Mpc) temperature map 
of Coma has been presented by Honda et al. (1996), who 
used GIS data from fourteen different ASCA pointings. 
Their map has a resolution of 30' (a GIS field of view) and 
reveals a hot region 40' northwest of the cluster as well 
as several other significant deviations from isothermality. 
Watanabe et al. (1997) presented a preliminary analysis of 
the temperature structure at smaller scales. 

In this paper, we use the ASCA observation of the Coma 
center to obtain a more detailed, ~ 5' resolution, temper- 
ature map of the central (r < 18'~ 0.7/i^q 1 Mpc) region 
of the cluster. All distance dependent quantities have as- 
sumed H = 50 km s _1 Mpc -1 , q = 0.5, all coordinates 
are given in the J2000 system, and unless otherwise noted 
all error bars are 90% confidence level. 

2. OBSERVATIONS & RESULTS 

ASCA observed the core of the Coma cluster on June 14, 
1993 (sequence # 80016000), and had useful exposures of 
7.6, 7.6, 9.4 and 9.4 ksec for the SIS 0-1 and GIS 2-3 detec- 
tors respectively. We have not included the thirteen offset 
pointings centered on this cluster (see Honda et al. 1996) 
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because they are subject to "stray light" contamination by 
photons scattered from the central area which we do not 
correct for in our analyses. 

Correctly characterizing the temperature distribution 
from ASCA data for an extended source on scales smaller 
than the entire field of one of the detectors requires a 
correction for the energy dependent PSF of the tele- 
scope (Takahashi et al. 1995). To this end, we have em- 
ployed two independent methods detailed elsewhere (Chu- 
razov et al. 1996, 1998, hereafter Method A; Markevitch 
et al. 1996a, 1998, hereafter Method B) to generate tem- 
perature maps with resolutions of ~ 5'. Method A pro- 
vides a continuous temperature map, although the correc- 
tion for the PSF is approximate (see Churazov et al. 1998 
for details). Method B derives spectra in projected regions 
defined by the user and takes exact account of the ASCA 
PSF. Note that both methods use the same ASCA PSF 
model and therefore may be subject to the same system- 
atic uncertainties. 

Figure 0a, developed using Method A applied to the 
ASCA GIS data, shows the approximate best fit tempera- 
ture, excluding the outer regions of the field of view where 
the uncertainty in the temperature was large. The map 
has been smoothed to a spatial resolution of ~ 5' to im- 
prove the signal-to-noise ratio. Contours in both panels of 
Figure show the filamentary structure in the X-ray inten- 
sity revealed by V97 after the removal of large-scale cluster 
emission via wavelet smoothing of the i?OS^4TPSPC im- 
age (see Vikhlinin et al. 1994 and V97 for details). In 
Figure 0b the wavelet contours are overlaid on the optical 
image from the Digital Sky Survey. We have added a box 
to delineate the area shown in Figure 0a. The positions 
of the two especially prominent galaxies, NGC 4874 and 
NGC 4889, are indicated in Figure 0a with solid trian- 
gles. There are three prominent features in the temper- 
ature map: two areas ~ 20% cooler than the average on 
the south and southeast edges of the map coincident with 
the filament, and a hot spot, ~ 25% hotter, displaced ~ 5' 
to the north from the peak of the intensity located ~ 1.5' 
southeast of NGC 4874. 

To better understand the statistical significance of the 
thermal deviations found in Figure 0, we employed Method 
B to produce a binned map of the temperature (Figure |2|) 
using data from both the GIS and SIS detectors. Integral 
to this method is the use of an intensity model which in 
this case we have chosen to be the ROSATPSPC image. A 
30" relative offset between the ROSAT and ASCA images 
was corrected so that model and data would be correctly 
aligned. The lowest wavelet contour from Figure has 
been overlaid in the figure in red, the ROSAT PSPC con- 
tours are overlaid in black, and the positions of NGC 4874 
and NGC 4889 are marked with filled triangles. The area 
displayed in Figure 0a is outlined for comparison. A plot 
of the temperature versus region is provided in the upper 
right panel of Figure 0. 

The integrated spectrum from within a 9' radius of the 
cluster center yields a best fit average temperature of 
9.0 ± 0.6 keV, consistent with previous results (Hughes 
et al. 1993, Honda et al. 1996). Figure shows that, 
except for regions 1, 9, 11 and 20, there are no signif- 
icant deviations from this average. Region 11 (kT = 
I2.7I2 Q keV) coincides with the center of the hot spot 



found in Figure 0, while regions 1, 9 and 20 (kT = 
3.91^0, 5. 7+^2 an d 6.lt^4 respectively) generally coin- 
cide with the cool regions along the filament in Figure 
As with previous work (Donnelly et al. 1998), the tempera- 
tures determined by both methods are in good agreement. 

3. DISCUSSION 

The simplest interpretation of our temperature maps, in 
agreement with B96, is that Coma is experiencing a merger 
with a moderately sized group associated with NGC 4874 
located south of the observed hot spot. The hot spot of 
emission is due to a bow shock generated by the group as 
it moves through Coma's intercluster medium(ICM). As 
the group has fallen toward the core of Coma, it has left 
a trail of cool gas delineating its in-bound track, which 
we see as the filament detected by V97. This is similar, 
for example, to simulations by Roettiger, Loken & Burns 
(1997, hereafter RLB97), Roettiger, Burns & Loken (1996, 
hereafter RBL96) and Norman & Bryan (1998, hereafter 
NB98). 

Recently both CD96 and B96 reported significant evi- 
dence for substructure in the core of Coma associated with 
the galaxies surrounding both NGC 4889 and NGC 4874. 
CD96 argued that the group around NGC 4874 was the 
"true central peak" , and the outer halo of the cD~ which 
would have been destroyed during core passage- has been 
re-aquired from merger debris during its long term resi- 
dence in the core. They were, however, hard pressed to 
explain the large velocity offset (~ 350 km/s) relative to 
the overall field required from their model for NGC 4874. 

Alternatively, the presence of a cD galaxy in the form 
of NGC 4874 could indicate that the group is in an early 
stage of the merger prior to entering the core, and as such 
has not yet been disrupted by the merger process. B96 sep- 
arated the cluster galaxies into "bright" and "faint" sam- 
ples, and found- within the "bright" sample- two galaxy 
groups, one associated with each of the two dominant 
galaxies. This is similar to the results of Fitchett & Web- 
ster (1987). B96 then argued that the "faint" sample was 
a better representative of Coma's overall velocity field and 
as such indicated that the two groups were "still in the 
process of merging with the cluster" . This model would 
explain the presence of the cD, the velocity offset for the 
NGC 4874 group as well as the temperature structures 
that we have found and the offset of the peak of Coma's 
overall emission from the two dominant galaxies (White 
et al. 1993, Vikhlinin et al. 1994). 

Following Markevitch et al. (1996b), we used the tem- 
perature maps to estimate the local relative specific en- 
tropy (entropy per particle) to better understand the gas 
merger process and the energetics involved. The relative 
specific entropy is defined as 



As = s — so = — fcln 
U 2 



T 



(1) 



where poi an< i s o are i m this case, the central values 
from an isothermal /3-model (p mo dei °c po(l + r 2 /r 2 )~ 3/3 / 2 , 
Cavaliere & Fusco-Femiano 1976). The emission- weighted 
average density in each projected region from the sym- 
metric /3-model was crudely corrected for the cluster's ob- 
served asymmetry by multiplying by the square root of 
the ratio of the actual brightness, as determined from the 
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Fig. 1. — Panel a shows the continuous temperature map produced using Method A (the scale bar is in keV). The map excludes regions 
where the uncertainty in the temperature is large, T + o~t /T — o~t > 1-5). Contours in both panels show the small scale structure in the Coma 
ROSAT PSPC image, detected by V97 after the removal of the main cluster emission using wavelet smoothing. In Panel b (reproduced from 
V97), the contours are overlaid on the Digital Sky Survey image. A square outlines the area from Panel a. The two central galaxies, NGC 
4874 and 4889, are on the right and left respectively and have been marked as triangles in Panel a. 




Fig. 2. — The binned temperature map produced by -the Method B. The black contours show the intensity from the ROSAT PSPC data and 
the red contour is the lowest contour shown in Figure hi Different colors indicate significant differences in temperature. The temperatures for 
each region with 90% uncertainties are plotted in the upper right panel. The lower right panel shows the estimated relative specific entropies 
for each region (see text for details). In the temperature panel, the average cluster temperature is shown as a dashed line. Dashed lines in the 
specific entropy panel show the entropy values for a symmetric, isothermal /3-model with T = 9 keV, relative to its central value, while the 
solid lines indicate average observed values in the groups, excluding regions 1 and 11 which are obviously discrepant. The positions of NGC 
4874 and NGC 4889, from right to left in the cluster core, are shown as filled triangles on the map, and a box has been provided to show the 
area displayed in Figure hk. 



ROSAT data, to the predicted model brightness: where p m odei is the value from the /3-model. This allows 

! us to compare the specific entropy for regions of interest 

/ ^obs \ 2 ,„% via the common benchmark of the /3-model, while taking 

P = Pmodei x (2) 

\ Pmodei J 



4 



A HOT SPOT IN COMA 



into account the smaller scale non-uniformities. We used 
the best fit /3-model parameters (3 = 0.75 and core radius 
r c = 10.5' from Briel et al. (1992). 

The lower right panel in Figure || shows the relative spe- 
cific entropy in each region with its 90% confidence error 
bar. The errors include contributions from the uncertain- 
ties in the surface brightness correction and the /3-model, 
but are strongly dominated by the contribution from the 
temperature determination. For clarity, we have segre- 
gated the regions into four groups, such that the relative 
specific entropy for the azimuthally symmetric isothermal 
model (shown as dashed lines) is constant within each 
group. The average relative specific entropy for each group 
(solid lines) is shown for comparison. Due to the obviously 
discrepant nature of both regions 1 and 11, they were ex- 
cluded from the average for their respective groups. 

In the core, Region 11 (the hot spot) has a significantly 
higher relative specific entropy compared to the other cen- 
tral regions. Region 8 also has a high relative specific en- 
tropy relative to its group, although this deviation is only 
marginally significant. Assuming that this hot gas lies in- 
side the core, and that initially it had a relative specific 
entropy equal to that of the rest of the core, we can esti- 
mate the total energy that has been injected into the gas 
as: 

Q ~ TAS ~ 0.3 - 1.0 X I0 el h~ hrgs. (3) 

By comparison, a gas mass of 2 x 10 12 Mq associated with a 
group falling from 5 Mpc into the core will acquire approx- 
imately 10-30 times this amount in kinetic energy, part of 
which it should lose prior to reaching the core. The con- 
version of the gas's gravitational potential energy is then 
a reasonable source for the energy required to produce the 
hot spot. 

We find that the temperatures in the cool regions (1, 
9 and 20) are consistent with the allowed lower limit of 
4 keV calculated by V97 for the filament. Regions 1 and 
(less significantly) 9 both also have low relative specific 
entropies compared to their group averages. Region 1 is 
especially notable in that it has a much lower relative spe- 
cific entropy than even the core values of the /3-model. The 
low values of the relative specific entropy for regions 1 and 
9 suggest that the gas in these regions only recently began 
to fall into the cluster and has not yet been completely 
mixed with or heated by the cluster's ICM. These regions 
are all spatially associated with the filament seen in X-ray 
image. 

In their analysis of the filamentary emission, V97 offered 
two possible explanations for the filamentary enhancement 
in emission. Either the filament arises from gas stripped 
via ram pressure from a group passing through Coma, or 
from a local perturbation in the ICM due to dark matter 
that has been tidally stripped from a merging group. The 
distinctive relative specific entropy of gas in the filament 
supports the former explanation. 

RLB97's simulations of a merger with a small group 
(see also RBL96 and NB98), find that the passage of a 
group through the cluster environment will lead to sub- 
stantial vorticity along the path of the in-bound object. 
This would temporarily introduce an additional pressure 
term in the in-falling gas preventing it from thermally equi- 
librating with the overall ICM. It is also possible that the 
vorticity would retard the mixing of the in-falling gas into 



the overall ICM, thus prolonging the lifetime of the fila- 
ment and addressing some the time scale issues raised in 
V97. 

On smaller scales, numerical simulations (RBL96, 
RLB97) indicate that, in a minor merger, the in-falling ob- 
ject will develop a protective bow shock, which will effec- 
tively shelter the galaxies from being ram pressure stripped 
of their gas. However, close to core passage the in-falling 
object's ICM encounters a rapidly increasing density in the 
gas of the main cluster. This would lead to the group's gas 
and galaxies decoupling from each other, and the galax- 
ies experiencing an "impulsive burst of ram pressure" as 
they pass through the shock. This would naturally lead 
to the galaxies being stripped of their gas prior to their 
coalescence with the cluster. Contrary to this, Vikhlinin 
et al. (1994) find a large gas halo associated with NGC 
4874, while Schombert (1988) finds the optical completely 
normal and unperturbed for a cD galaxy. 

Radio observations of the galaxy NGC 4869, ~ 4' south- 
west of NGC 4874, by Feretti et al. (1990) indicate that 
it is in orbit around the cD galaxy. They noted a ma- 
jor problem in the confinement of the radio tail associated 
with NGC 4869 due to its internal pressure strongly de- 
creasing with distance from the center of the galaxy. One 
possible explanation is that, due to a projection effect, 
the pressure of the ambient medium is lower than that 
which they have computed. Feretti et al. (1990) reject 
this based upon the assumption that NGC 4874 is nearly 
at rest at the bottom of the cluster potential. If, however, 
NGC 4874 and NGC 4869 are part of a group still falling 
into Coma a significant projection effect may be at work. 
Given the large extent of the X-ray gas- more than 12' in 
diameter- found by V97 around NGC 4874, it is plausible 
that NGC 4869 is being shielded from the overall cluster 
ICM (RLB97), and the relevant ambient medium is that 
of the in-falling group. 

These arguments suggest that the group has not quite 
arrived at the core. Further, the velocity studies of Mellier 
et al. , B96 and CD96 find a statistically significant offset 
for the NGC 4874 group relative to the Coma's average 
velocity field. This suggests that that the group is falling 
in from the front side of Coma at a mild inclination to our 
line of sight. 

There are some complications to this simple model due 
to the nearby presence of another subcluster associated 
with NGC 4889, ~ 7' to the east from NGC 4874. Baier 
et al. (1990) find that the group of galaxies associated 
with NGC 4889 is elongated along an axis from the south- 
east to the northwest, possibly indicating a recent dynam- 
ical/tidal interaction similar to that described by RLB97 
as occurring after core passage. It is possible that the 
slightly high entropy in region 8 briefly mentioned above 
is related to this activity. The lack of an extended envelope 
around NGC 4889 also would be expected if the group has 
recently passed through the core. As with the NGC 4874 
group, CD96 and B96 find a statistically significant veloc- 
ity offset for the NGC 4889. This suggests that both the 
NGC 4889 and NGC 4874 groups - in that order- may be 
recent arrivals in the core. We are thus seeing the debris 
of the final stages of one minor merger, just as another 
minor merger is occurring. 

Other complications arise from what the model de- 
scribed above suggests that we should see. Nominally, we 
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expect that the surface distribution of galaxies in an in- 
falling group would experience some distortion prior to 
core passage. In their simulations, RLB97 show just such 
a morphology at -0.25 Gyr for a 2:1 mass ratio merger, 
although the distortion is less for the secondary (i.e. in- 
falling) cluster than for the primary and it is not clear 
how this effect would scale with mass. However, Baier 
et al. (1990) find that the group of galaxies around NGC 
4874 is circularly distributed suggesting no dynamical per- 
turbations. This lack of distortion evident in the NGC 
4874 group may be due to its being still quite far from 
core passage. 

Finally, we would expect the temperature and entropy 
results for Region 5- situated directly between regions 1, 
9 and 20- to share in the properties of the filament as out- 
lined above. However, we find that it is nearly identical to 
the non-filament gas. Perhaps this is due to a projection 
or dilution effect due to the increasing brightness of the 
overall ICM. Unfortunately, the resolution of our data is 
inadequate for further speculation. 

4. SUMMARY 

We have analyzed the ASCA pointing covering the cen- 
tral r=18' area of the Coma cluster. Although, the tem- 
perature map of the gas in the core of the cluster does 
not exhibit the strong temperature variations observed in 
clusters undergoing major mergers, we find several signif- 
icant features in the temperature map indicating recent 
and ongoing mergers with small groups. 

Our analysis indicates that gas to the south/southeast 
side of the cluster is significantly cooler and has a lower 
relative specific entropy than the average. The regions 
of low temperature and relative specific entropy coincide 



with the gas filament, detected by V97, extending from 
the cluster's center to the south and east towards NGC 
4911. This low entropy gas is most likely deposited in the 
wake of the NGC 4874 group as it falls toward the cluster. 

The temperature map also reveals a hot spot located to 
the north of the cD galaxy, NGC 4874, and the peak of 
the intensity profile. This is most likely a bow shock in the 
gas which has developed as this small group of galaxies has 
fallen in to the core of the cluster. Our rough estimate of 
the relative specific entropy of the shocked gas is consis- 
tent with heating from the conversion of the kinetic energy 
of a small group falling into the cluster. 

The X-ray and optical data suggest that a small group 
of galaxies associated with NGC 4874 is in-bound towards 
the cluster core from the southeast leaving a tail of cooler 
gas along its trajectory. The relatively unperturbed state 
of the group suggests that core passage has not yet oc- 
curred. The group associated with NGC 4889, although 
still in the process of being incorporated into the cluster, 
has already passed the core and most likely been resident 
in the core for some time, as evidenced by its dynamically 
disturbed state. 

Future observations of Coma with ^IXAFwill allow us to 
study the characteristics of the ICM around these merger 
remnants at much higher spatial resolution. This will help 
us to understand the processes at work during and after 
mergers in much greater detail. 

We would like to thank the referee for their helpful sug- 
gestions and comments. RHD, MM, WF, CJ acknowledge 
support from the Smithsonian Institute and NASA con- 
tract NAS8-39073. 
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